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Abstract. We discuss the problem of counting incidence matrices, i.e. zero-one matrices with 
no zero rows or columns. Using different approaches we give three different proofs for the leading 
asymptotics for the number of matrices with n ones asn-> oo. We also give refined results for 
the asymptotic number of i x j incidence matrices with n ones. 



1. Introduction 

We call an incidence matrix a zero-one matrix with no zero rows and columns and denote by 
F(ri) the number of incidence matrices with exactly n ones, where n £ N. For example, the four 
incidence matrices with n = 2 are 



taken from the On-Line Encyclopedia of Integer Sequences [7], where this appears as sequence 
A101370. For convenience, we further define F(0) = 1. 

If one imposes additional symmetries or constraints, such as allowing or prohibiting repeated 
rows or columns, or considering equivalence classes under row or column permutations, one is 
led to many different enumeration problems, as discussed in jSJ. 

The counting problem can be interpreted in a surprisingly rich variety of different ways, 
leading to rather different mathematical approaches. 

• Counting hypergraphs by weight 

Given a hypergraph on the vertex set {x±, . . . , x r }, with edges E\, . . . , E s (each a non-empty 
set of vertices), the incidence matrix A = (djj) is the matrix with entry 1 if Xi G Ej, 
and otherwise. The weight of the hypergraph is the sum of the cardinalities of the edges. 
Thus F(n) is the number of vertex- and edge-labelled hypergraphs of weight n with no 
isolated vertices, up to isomorphism. 

• Counting bipartite graphs by edges 

Given a zero-one matrix A = (Aij), there is a (simple) bipartite graph whose vertices are 
indexed by the rows and columns of A, with an edge from rj to Cj if A^ = 1. The graph 
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The first few terms of the sequence F(n) for n £ N are 



1, 4, 24, 196, 2016, 24976, 361792, 5997872, 111969552, . . . 



has a distinguished bipartite block (consisting of the rows). Thus, F{n) counts labelled 
bipartite graphs with n edges and a distinguished bipartite block. 

• Counting pairs of partitions, or binary block designs 

A block design is a set of plots carrying two partitions, the treatment partition and the 
block partition. It is said to be binary if no two distinct points lie in the same part of both 
partitions; that is, if the meet of the two partitions is the partition into singletons. Thus, 
F(n) is the number of binary block designs with n plots and labelled treatments and blocks. 

• Counting orbits of certain permutation groups 

A permutation group G on a set X is oligomorphic if the number F*{G) of orbits of G on 
re-tuples of elements of X is finite for all n. Equivalently, the number F n (G) of orbits on 
ordered re-tuples of distinct elements of X is finite, and the number f n (G) of orbits on re- 
element subsets of X is finite, for all re. These numbers satisfy various conditions, including 
the following: 

• K(G) = ELi S ( n > k ) p k{G) and its inverse F n (G) = £™ =1 s(n, k)F*(G), with s(n, k) 
and S(n, k) Stirling numbers of the first and second kind, respectively; 

• f n {G) < F n (G) < n\f n {G), where the right-hand bound is attained if and only if the 
group induced on a finite set by its setwise stabiliser is trivial. 

For example, let A be the group of all order-preserving permutations of the rational numbers. 
Then f n (A) = 1 and F n (A) = n\ . 

Now if H and K are permutation groups on sets X and Y, then the direct product 
H x K acts coordinatewise on the Cartesian product X x Y. It is easy to see that 



Let [x\, yi), . . . , (x n , y n ) be re distinct elements of X x Y . If both X and Y are ordered, 
then the set of re pairs can be described by a matrix with re ones in these positions, where 
the rows and columns of the matrix are indexed by the sets {x%, . . . , x n } and {yi, . . . , y n } 
respectively (in the appropriate order). Thus 



Discussion of this "product action" can be found in [3]. 

For an extended discussion of these interpretations see [JJ]. For instance, when considering 
hypergraphs it is more natural to consider the unlabelled problem, which leads to identification of 
incidence matrices which are equivalent under permutation of rows or columns. Also, forbidding 
repeated rows corresponds to counting simple hypergraphs with no repeated edges. 

2. The asymptotics of F(n) 

It is possible to compute F(n) explicitly. For fixed re, let my(n) be the number of i x j matrices 
with n ones (and no zero rows or columns). We set moo(0) = 1 and F(0) = 1. Then 



F*(H x K) = F*(H)F*(K). 



F(n) = f n {A x A). 




(1) 



so by Mobius inversion, 




and then 




(3) 



i<n j<n 



For sequence a n , b n , we use the notation a n ~ b n to mean lim n ^oo a n /b n = 1. It is clear from 
the argument above that 

F{n) < ( n2 \ ~ -^=(ne)», 
and of course considering permutation matrices shows that 

Fin) > n\ ~ \j2irn ( — 
V e 



Theorem 1 
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We remark that for n = 10, the asymptotic expression is about 2.5% less than the actual 
value of 2324081728. 

As announced in [5], we have three different proofs of Theorem^ The first proof employs 
pairs of random preorders and a probabilistic argument, the second proof uses counting of orbits 
of products of permutation groups, and the third proof employs a surprisingly simple identity. 
First Proof: This proof uses a procedure which, when successful, generates an incidence matrix 
uniformly at random from all incidence matrices. The probability of success can be estimated 
and the asymptotic formula for F(n) results. 

Let R be a binary relation on a set X. We say R is reflexive if (x, x) £ R for all x £ X. We 
say R is transitive if (x, y) £ R and (y, z) £ R implies (x, z) £ R. A partial preorder is a relation 
R on X which is reflexive and transitive. A relation R is said to satisfy trichotomy if, for any 
x, y £ X, one of the cases (x, y) £ R, x = y, or (y, x) £ R holds. We say that R is a preorder if 
it is a partial preorder that satisfies trichotomy. The members of X are said to be the elements 
of the preorder. 

A relation R is antisymmetric if, whenever (x, y) £ R and (y, x) £ R both hold, then x = y. 
A relation R on X is a partial order if it is reflexive, transitive, and antisymmetric. A relation 
is a total order, if it is a partial order which satisfies trichotomy. Given a partial preorder R 
on X, define a new relation S on X by the rule that (x, y) E S if and only if both (x, y) and 
(y, x) belong to R. Then S is an equivalence relation. Moreover, R induces a partial order x on 
the set of equivalence classes of S in a natural way: if (x,y) £ R, then (x, y) £ R, where x is 
the S'-equivalence class containing x and similarly for y. We will call an S-equivalence class a 
block. If R is a preorder, then the relation R on the equivalence classes of S is a total order. See 
Section 3.8 and question 19 of Section 3.13 in for more on the above definitions and results. 
Random preorders are considered in jH]. 

Given a preorder on elements [n] := {1,2,..., n} with K blocks, let B±,B2, . . . , Bk denote 
the blocks of the preorder. Generate two preorders uniformly at random, B±, B2, ■ ■ ■ , Bk and 
B[, B' 2 , . . . , B' L . For each 1 < i < j < n, define the event Dij to be 

Dij = {for each of the two preorders i and j are in the same block}. 

Furthermore, define 



W= J A 

l<i<j<n 

where the indicator random variables are defined by 

Id 



1 if occurs; 



ij I otherwise. 



If W = 0, then the procedure is successful, in which case B k n B[ consists of either or 1 
elements for each 1 < k < K and 1 < I < L. If the procedure is successful, then we define the 
corresponding K x L incidence matrix A by 



1 if B k n B[ £ 0; 
if B k n B{ = 0. 



It is easy to check that the above definition of A in fact produces an incidence matrix and that 
each incidence matrix occurs in n\ different ways by the construction. It follows that 

_ P(n?nW = 0) 
F{n) ~ n\ ' 

where P(n) is the number of preorders on n elements if n > 1 and P(0) = 1. 

It is known (see for example) that the exponential generating function of P(n) is 

f^ n = — • (4) 
^ n\ 2-e z y ' 

n=0 

The preceding equality implies that P{n) has asymptotics given by 

P(n) ~ ^ ( r^)^ 1 ■ (5) 



2 V>g2 

It remains to find the asymptotics of P(W = 0). 
The rth falling moment of W is 

E(W) r = EW(W- 1)---(W -r + 1) 

= El Y, l '.:J:---l' :J ) ( 6 ) 

\ pairs (is,j s ) different 
\ aii i s and j s different 

* 

with ^2 defined to be the sum with all pairs (i s ,js) different, but not all i s ,j s different. 

First we find the asymptotics of the first term in (J7J). For given sequences i%, 12, ■ ■ ■ , i r , 
ji,j2, ■ ■ ■ ,jr, the expectation ■ ■ ■ h rt j r ) is the number of ways of forming two preorders 

on the set of elements [n] \ {ji, J2, • • • , jr} and then for each s adding the element j s to the block 
containing i s in both preorders (which ensures that Di s j g occurs for each s) and dividing the 
result by P(n) 2 . Since the number of ways of choosing i\ , 12, ■ ■ ■ , i r , ji, ji, ■ ■ ■ , j r equals 2 r (n-2r)\ > 
This gives 



, aii i a and j s different 



n\ P(n — r) z 
2 r (n - 2r)\ P(n) 2 

(log2)^ r 



where we have used © • 



The second term is bounded in the following way. For each sequence {12^32)1 ■ ■ ■ , (is,js) 

in the second term we form the graph G on vertices (Js=lOs> 3*} W1 th edges Us=i{{*s> is}}- Con- 
sider the unlabelled graph G' corresponding to G consisting of v vertices and c components. The 
number of ways of labelling G' to form G is bounded by n v . The number of preorders corre- 
sponding to this labelling is P{n — v + c) because we form a preorder on n — v + c vertices after 
which the vertices in the connected component of G containing a particular vertex get added to 
that block. Therefore, we have 




P{nf 
£o(n 2 -), 



G> 



where the constant in O (n 2c ~ v ^ is uniform over all G' because v < 2r. Since at least one vertex 
is adjacent to more than one edge, the graph G is not a perfect matching. Furthermore, each 
component of G contains at least two vertices. It follows that 2c < v and, as a result, 



E \J2 I iuh--- I ir,jr) =0{n 

The preceding analysis shows that 

E(W) r r 



— 1> 



2\ r 



(log 2) 
2 

for each r > 0. The method of moments implies that the distribution converges weakly to the 
distribution of a Poisson((log 2) 2 /2) distributed random variable and therefore 

F(W = 0) ~ exp M lQ g 2 H . ( 8 ) 



□ 

Second Proof: We now give a proof using product actions of groups, as discussed in the 
introduction. First of all, this approach leads to a different and simpler expression than © for 
F(n) as a sum of terms of alternating sign. 

Proposition 2 



n 

M.I ^— ' 

where 



k=l 



P(n) = ^2S{n,k)k\ 

k=l 

is the number of (total) preorders of {1, . . . , n}, and s(n, k) and S(n, k) are Stirling numbers of 
the first and second kind respectively. 

This is proved in [2], but can be seen as follows. Using the group A of all order-preserving 
permutation groups acting on Q, we consider the direct product A x A acting on Q x Q. We 
have F n (A) = n\, whence it follows that F*(A) = Y2=i S ( n i k ) kl = P ( n )- Thus 



n 

F*{A xA) = P(n) 2 = Y, S(n, k)F k {A x A), 

k=l 



and so the inverse relation between the two kinds of Stirling numbers gives 



n 

F n (A x A) = s(n, k)P(k) 2 . 

k=l 

Finally, the group A x A has the property that the setwise stabiliser of a finite set fixes it 
pointwise, and so f n (A x A) = F n (A x A)/n\ .□ 

We now replace P(k) by the asymptotic form (J5J) given earlier. For k > n/2, the difference is 
exponentially small; and we will show below that the contribution of the terms with k < n/2 is 
negligible, so it suffices to note that the error we make is smaller than the approximated term. 

So let 

1 1 n 

F'(n) = -.-Y,s(n,k)(k\fc^\ 
' k=l 

where c = l/(log2) 2 is as in the statement of the theorem. As we have argued, F(n) ~ F'(n). 

Now (— l) n ~ k s(n, k) is the number of permutations in the symmetric group S n which have k 
cycles. So we can write the formula for F'{n) as a sum over S n , where the term corresponding to 
a permutation with k cycles is (— l) n ~ k (k\) 2 c k+1 . In particular, the identity permutation gives 
us a contribution 

g(n) = \n\c n+ \ 

and we have to show that F'(n) ~ Cg(n) as n — > oo, where C = exp(— (log 2) 2 /2). 

To prove this, we write F'(n) = F[(n) + F^n) + F^(n), where the three terms are sums over 
the following permutations: 

F[: all involutions (permutations with a 2 = 1); 
i^: the remaining permutations with k > [n/2]; 
Fg: the rest of S n . 
We argue that F[(n) ~ Cg{n), while F^n), F^n) = o{g{n)). 



Case F[: Let I = n — k. Now an involution with k cycles has I cycles of length 2 and n — 21 
fixed points; so I < n/2. The number of such permutations is 



So 



n\ (2l)\ _ n{n - 1) ■ • • (n - 21 + 1) 
2l)Yn ~ ¥T\ 



F'M ^ n {n-l)...{n-2l + l) y ((n-/)!) 2 , 

g(n) ^ 2H\ K ] (n!) 2 

L " /2J 1 f-l\ l (n-l)---(n-2l + l) 



E 



Now 



, l\\2cj n---(n-l+l) 



^ 11 V 2c J 6XP V 2c 



i=0 

so we have to show that the factor involving n makes no difference to the limit. Now this factor 
is always less than 1, so the series is absolutely convergent (and uniformly in n); so we can 
choose r large enough that the sum of r terms of each sequence is close to its limit. Then, since 
the factors tend to 1 as n — > oo, for n large each of these r terms is close to its limit. So the 
assertion is true: that is, F[(n) ~ Cg(n). 



Case F^: A permutation which has k = n — I cycles and is not an involution has at least 
n — 21 + 1 fixed points, and there are at most 



n 

21 - 1 



(21 - 1)1 = n{n - 1) • • • (n - 21 + 2) 



such permutations. So, ignoring signs, 

F 2( n ) sr^{n-l){n-l-l)---(n-2l + 2) _ t 
g{n) ~ n(n-l)---(n-l + l) 

1 1 
< - 



n 1 — c 1 
which is 0(l/n). 

Case F%: We simply observe that there are at most n! such permutations, so 

\n/2\ 



gin) 



k=i v 7 



Now n\/(k\) 2 > (i n y 2 j) — — e )™ ^ or l ar § e n > so this sum is 0((2 — e) n n/2) = o(l) as n — » oo. 
□ 

Third Proof: If one is interested in asymptotic enumeration of F(n), (J2J), being a double sum 
over terms of alternating sign, is on first sight rather unsuitable for an asymptotic analysis. The 
expression in Proposition |2] is also an alternating sum. We present a derivation of the asymptotic 
form of F(n) based on the following elegant and elementary identity, which gives F(n) as a sum 
of positive terms. 



Proposition 3 



Proof. Insert 



k=0 1=0 



(9) 



k=i v 7 l=j KJ/ 



into ((SJ) and resum using Q.D 

We start the asymptotic analysis by rewriting © as 

»^») = EE^f|. m 

where (x) n = x(x — 1) . . . (x — n + 1) is the falling factorial. Given the identity 

P(n) = £ ^TT < ( 12 ) 

fc=0 

which follows from expanding (J3J, (|TTj) is bounded above by P(n) 2 , as the factor {kl) n / (kl) n 
takes values in [0, 1]. 



For n < kl, a straightforward expansion of the factor gives 



(fcZ)n 
(kl) n 



exp 



™ B j+1 (n) - B j+1 (0) 



(13) 



Here, we have used that Ylk=o ^ = (-^i+i( n ) ~~ -Bj+i(0))/(j + 1) where Bj{x) is a Bernoulli 
polynomial. It follows that 



{|^ = e-£ (l + 0(n/W) + 0(n 3 /(W) a )) ■ 



(14) 



(This argument will be presented more thoroughly for (z) n with complex-valued z in the next 
section.) The sum (jllj) is dominated by terms around k = I = n/log2, so that we expect 
the correction to give e~( log2 ) / 2 , which in turn would imply n\F(n) ~ P(n) 2 e _( - log2 ) / 2 . The 
difference is given by 



!F(n) - P(n) 2 e-^ 2 ) 2 / 2 = -^Jl (gk _ e -(.o g2 )V 2 

fc=0 z=o \\ J 



(15) 



To proceed we choose mo < n/log2 < mi and split the summation. We obtain 



\F{n)-P{nfe-^ 2 ' 2 



7711 mi /c n P 

£ E E 

k=rriQ l=m,Q 



(kl) 



(kl) 



n _ e _(log 2) 2 /2 



+ 2 E 



A:'' 



2 fc+i 



mo — 1 x 



fc=0 

Specifying mo = n./ log 2 — cn s and mi = nj log 2 + en* 5 for 1/2 < 5 < 1 and c > 0, we use (|14|l 
to estimate 



k n 

2 k + 1 1 2 fc+1 

fc=0 fc=mi + l 



(16) 



« e -(log2) 2 /2 



e-^i (1 + 0(n/kl) + 0(n 3 /(kl) 2 )) - e -( lo s 2 ) 2 / 2 
e -(lo 6 2)V2 ^ + oCn*" 1 )) (1 + 0(n- x )) - e -( lo s 2 ) 2 /2 



0(n 



for mo < k, I < mi. This allows us to bound the first term in (|16|) by P(n) 2 0(n s 1 ). To get a 
bound on the second term, we utilize the following 

Lemma 4 (a) For K,n £ N and K < n/ log 2, 



A A: n K n e n ' K 1 



fe=0 



2 fc+i - 2^ 2 e n l K - 2 ' 



(17) 



(b) For K,n G N and if > n/log2, 



« fc n K n e n/K j 
2^ 9^+1 



2 fc+i - 2^ 2 2 - e"/-^ ' 



(18) 



Proof Part (a) follows from the estimate 

5^2^+1 = 2^+1 ^ ~2 k ~ Z]r ~ ~ X) 2 *^ 1 ~~ k / K ) n ^ -^kTT ^2 \ 2e n K ) 

k=0 k=0 k=0 k=0 

and part (b) similarly from 

k n K n ^ {k/K) n < K n ^(l + k/K) n < K n ^fe n / K \ k 

k=K+l k=K+l k=l k=l \ / 

□ 

For K = nj log 2 =p cn s , we find 

r^n n/K i n -n , 

^£ 1 = ^-—e~ an ^Oin 1 - 5 ) = P(n)0 e 

2 K 2 \e n / K -2\ (log 2)™ v ; v ' V 

where a = c 2 (log2) 2 /2. Using Lemma we therefore bound the second term in (j!6ft by 
P(n) 2 (exp(— an 2 * -1 )) . Altogether we find 

n\F{n) - P(n) 2 e- (log 2)2/2 = P(n) 2 (o(ra 5_1 ) + O (expf-a/- 1 ) 

and as 1/2 < 6 < 1, we have 



n!-F(n 

mil 

n— >oc 

which completes the proof. □ 



n^nj _ ( i og 2) 2 /2 



3. The asymptotics of rriki{n) 

In this section we present results on the number of incidence matrices with specified numbers 
of rows and columns. To obtain the desired asymptotic form of m k i{n) from eqn. ((2j), we need 
to deal with the challenge that summing over large terms with alternating signs can lead to 
enormous cancellations. Fortunately, there is a standard trick using the calculus of residues. 

Proposition 5 

m k i{n) = — Res(-- — ; s = oo, t = oo) . (19) 

n! (s) k+ i[t)i +1 

Proof Using the fact that 

(-1)™ 

Res(r(s), s = -m) - 



we write 



m k i{n 



ml 



nl { > j^j^ i\ j\ (k-i)\(l-j)\ 

k w. (-i)^ f r (st )n 



n\ (2vri) 2 J Cl _ kQ] J C[ _ io] ( s + k) k+1 (t + l) l+1 
kU[ 1 f dsf dt 



n\ (2vri) 2 y C[o ^ 7 C[o i] (s) k+1 (t) l+1 



Here, the contours C[ a ,fc] encircle the (real) interval [a, b] counterclockwise. As the integrand is a 
rational function, the contour integrals can be expressed as residues at infinity. □ 



This formulation allows us to do the asymptotic analysis via saddle point analysis of a contour 
integral. We consider the scaling behaviour of mk,i(n) as n — > oo with k = nn and I = An for 
fixed A, k. As a preparation, we state the following Lemma. 



Lemma 6 Let n S N and z E C with \z\ > n. Then 

B j+1 (n) - B j+1 (0) 



z exp 



§ A/ • 



(20) 



Moreover, we have the asymptotic expansion 



log(,z) n ~ (z + 1/2) log z — (z — n + 1/2) log(z — n) — n 

+ ^2fc(2fc-l) (i2fc=I " ( z _ n )2fe-i) ' ( 21 ) 

as |z — n| and |z| tend to infinity. Here B n {x) is the n-th Bernoulli polynomial and B n = B n (0) 
the n-th Bernoulli number. 

Proof. We write 

oo 



w.-»-n(i-i) 

fc=0 V 7 



For \z\ > n we take logarithms and expand log(l — k/z) in k/z. Exchanging the order of 
summation and using that 



n-l 



^ = (B j+1 (n) - B j+1 (0))/(j + 1) 



k=0 

gives (|2()|). One can obtain (f2*T|) by substituting B n (x) = X^fe=o (k)BkX n ~ k and exchanging the 
order of summation again. As the double sum here is not absolutely convergent, the resulting 
series (|21|) cannot be expected to be convergent. Instead of labouring to prove that one still 
arrives at an asymptotic expansion, we point out that the result is just the difference between 
the Stirling series for log(z!) and log((z — n)!) and argue that standard arguments used in the 
derivation of the Stirling series also lead to (|21|) . In contrast with the Stirling series, which is 
valid for | arg(z)| < n, the validity of this expansion is not restricted to a sector of the complex 
plane. □ 

From (1201) we obtain that 



(z) n = z n e-t (1 + 0{n/\z\) + 0(n 3 /\z\ 2 )) , (22) 

whereas from (|21|) we obtain that 

(z) n = (z- nr-^/V+^e-™ (1 + 0(l/\z\) + 0(l/\z - n|)) . (23) 

We now state the main theorem of this section. 
Theorem 7 For fixed a, r > 0, 

rj 2n i 

m k i(n) -e nw{a) v(a)e nw{T) v(T)e-^ (24) 

re! 



with 

w(x) = x{l-e- l l x )\og{l-e- l l x ) + \ogx-e- l / x 
v{x) = 1/ 



x(l - e~ l / x ) 



x(l - e- l l x ) - e- x l x 

and k = na(l - e" 1/,J ), I = nr(l - e~ 1/r ). 
Proof In order to evaluate the integral 

nl ( \ 1 f A [ M ^ 

:m k i(n) = . ,. Q / ds dt 



kill ' (2vr*) 2 J C[0M Jc m Wfc+l(t) 



H+l 



asymptotically, we approximate the integrand uniformly using (|22[) and (|23j) on contours 
satisfying \s\ = R s > k and \t\ = Rt > I- We find 

n! , , e fc +' f , /■ (a - A;)'-*- 1 ^ (t - J)*-I-V2 „ 2 

if" 1 "' = W J\s\=r s ds L- R , * s-'^n ^ + *> 

where 

i? = 0(n/i? SJ R t ) + 0(n 3 /{R s R t ) 2 ) + 0(l/i? s ) + 0(l/i? t ) + 0(l/|i2 a - fc|) + 0(l/|i^ - Z|) . 
Substituting s = ncr, t = nr, k = nK, I = nX, this simplifies to R = 0(l/n) and we arrive at 



m kl (n) ^ K nK X nX n n e n ^^ I da I dr e nf( - a ^g(a, K )e nf( - T ^ g(r, A)e"2^ 

(2my y M=P(T j 1t{=Pt 

with R s = np a and Rt = np T and 

f(x,y) = (x - y)log(x - y) - (x - l)logx and y) = ((x - y)x) _1/2 . 

As n tends to infinity, each integration is dominated by saddle points <j s and t s on the positive 
real axis. The saddle point equation is 

= d 1 f(a s ,K) = ]og(l-—) +- , 

(here, d\ denotes taking the derivative with respect to the first argument) and an identical 
expression for t s . There exists a unique positive solution ct s (k), and a standard saddle-point 
evaluation (see e.g. [2]) gives 



m«(n) ~ fi «A"Ve^e"«^ K ) 5 (, SlK )e"/(--^(r Sl A)e-^- ^ 



which simplifies to the desired result. □ 

Theorem d can be used for a fourth proof of Theorem ^ via an asymptotic evaluation of the 
sum over niki (n) . The sum is dominated by terms near a s = t s = 1/ log 2 from whence it follows 
that the distribution has a peak about k s = l s = n/(21og2). 

We conclude this paper with giving an identity for m/y(n) which is a refinement of Proposition 

ca 



Proposition 8 




(25) 



Proof. 

Using (ij) n = Y^=\ s ( n > r )(^jYi we write (j2J) as 



m k i(n) = — ] Yl s ( n > r ) Yl 




and resum using k\S(r, k) 
Inversion of ()25|) gives 



£,< fc (-i) fc - l (tK- n 



n 



k\S(n,k)l\S(n,l) = ' s ^r\S{n,r)m] i i{r) 



(26) 



r=l 



which has a straightforward combinatorial interpretation, as r\S(n, r) is the number of preorders 
of n elements into k blocks. The left hand side of (|26|) is just the number of ways of choosing two 
preorders of an n-set into k and / blocks, respectively. The right hand side of (|26j) corresponds 
to counting the number of ways in which elements of an n-set can be distributed into r cells of 
a k x /-array, where the cells are given by k x /-incidence matrices with r ones, for arbitrary r. 

Summing (|25|) over k and / provides another proof of Proposition [21 Proposition El could also 
be used as a basis for Theorem We leave this as an exercise for the reader. 
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